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The neutron spin rotation induced by parity-violating (PV) components in the nucleon-nucleon 
potential is studied in n-d scattering at zero energy. Results are obtained corresponding to the 
Argonne vig two-nucleon and Urbana-IX three-nucleon strong-interaction potentials in combination 
with either the DDH or pionless EFT model for the weak-interaction potential. We find that this 
observable is dominated by the contribution of the long-range part of the PV potential associated 
qq ■ with pion exchange. Thus its measurement could provide a further constraint, complementary to 

' that coming from measurements of the photon asymmetry in n-p radiative capture, on the strength 

, of this component of the hadronic weak interaction. 
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I. INTRODUCTION 

m 

Cn| I We report on a calculation of the neutron spin rotation in n-d scattering at zero energy. The present work is the 
continuation of a program to study parity-violating (PV) effects in few-nucleon systems induced by hadronic weak 
interactions. Two earlier papers dealt with the two-nucleon systems: the first [l[ was devoted to p-p elastic scattering, 
and presented a calculation of the longitudinal asymmetry in the lab-energy range 0-350 MeV. The second [2| provided 
a detailed account of the i) neutron spin rotation in n-p scattering at zero energy, ii) longitudinal asymmetry in n-p 
elastic scattering up to lab energies of 350 MeV, iii) photon angular asymmetry in n-p radiative capture at thermal 
neutron energies, and iv) photon helicity dependence of the d(y, n)p cross section from threshold up to energies of 20 
MeV. 

In this study we adopt two different models for the PV potential. One is parameterized in terms of tt-, p-, and 
w-meson exchanges, in which the weak couplings of these mesons to the nucleon are estimated within a quark-model 
approach incorporating symmetry arguments and current algebra requirements. It is known as the DDH model [H, 
and has provided for over two decades a useful framework for the analysis and interpretation of measurements of PV 
, observables (see in this connection the recent review in Ref.[4j]). 
£T) ■ The other model is that derived recently by Zhu et al. [5j in an effective field theory (EFT) approach, in which 
\ only nucleons are retained as explicit dynamical degrees of freedom, while pions and baryon resonances, such as A 
isobars, are integrated out. It is formulated in terms of a number of four-nucleon contact terms. Obviously, its use is 
qq ' restricted to processes occurring at energies much below the pion mass. 

The rotation of the neutron spin in a plane transverse to the beam direction is calculated perturbatively to first order 
in the DDH and EFT PV potentials. We use n-d wave functions obtained with the hyperspherical-harmonics (HH) 
method @, 0, @ from realistic strong-interaction Hamiltonians, consisting of the Argonne v ls H (A V18) two-nucleon 
potential with and without the inclusion of a three-nucleon potential, the Urbana IX model [10] . The AV18/UIX 
Hamiltonian quantitatively and successfully accounts for a wide variety of few-nucleon bound-state properties and 
reactions, ranging from binding energies and charge radii to electromagnetic form factors, low-energy scattering 
observables and electroweak capture cross sections (for a review, see Refs. [TP. Il2j). 

Our prime objectives in this as in the earlier papers are to develop a systematic and consistent framework for 
studying PV observables in few-nucleon systems, where accurate microscopic calculations are feasible, and to use 
available and forthcoming experimental data on these observables to constrain the strengths of the long- and short- 
range parts of the two-nucleon weak interaction. While no measurements of the neutron spin rotation in n-d scattering 
have been carried out so far to the best of our knowledge, nevertheless we hope the present study will provide strong 
motivation for undertaking this type of experiments. Within the context of the calculations based on the DDH 
model, we find that this observable is dominated by the contribution associated with the long-range pion-exchange 
component. Indeed, assuming a liquid deuterium density of p = 0.4 x 10 23 atoms cm -3 and the "best" values of Ref. Q 
(as listed in Table II of Ref. [2|) for the tt-, p-, and cj-meson weak (and strong) coupling constants to the nucleon, the 
predicted spin rotation per unit-length of matter traversed is 0.52 x 10 -7 rad-cm -1 with the complete DDH potential, 
and 0.46 x 10" 7 rad -cm when only its pion-exchange term is retained. We find the model dependence of this result 
upon the (realistic) strong-interaction potentials used to generate the n-d continuum states to be quite small. Thus, 
a measurement of the n-d spin rotation could provide a further constraint, complementary to that coming from the 



2 



measurements in progress of the photon asymmetry in the n-p radiative capture, on the strength of the long-range 
part of the PV potential. 

The remainder of this paper is organized as follows. In Sec. [TT] the neutron-spin-rotation observable is related to 
matrix elements of the PV potential between incoming and outgoing n-d states. The DDH and EFT PV potential 
models are briefly discussed in Sec. IIIII in particular, in Appendix [Al it is shown that the number of independent 
four-nucleon contact terms in the EFT formulation is five rather than ten, as obtained in Ref. [5j. The Monte Carlo 
techniques used to calculate the relevant matrix elements are reviewed in Sec. IIV1 where the details of a test calculation 
based on simple Gaussian wave functions are also described. Finally, Sec. [V] contains a brief discussion of the results 
along with some concluding remarks. 



II. THE NEUTRON SPIN ROTATION OBSERVABLE 



In first order perturbation theory in the weak interactions, the neutron spin rotation per unit-length of matter 
traversed, d(j)/dz, in n-d scattering is given by pi. [l3T| 



"~T~ = T~~ — ^ c ( '(PZjmn, m d | v PV | pi;m n , m d ) (+) , (2.1) 

pdz 3 v IC \ ^ 



where p is the density of deuterons (a parameter under control of the experimenter), w PV denotes the PV nuclear 
potential, | pi] m n , ra^fr' are the n-d scattering states with incoming-wave (— ) and outgoing-wave (+) boundary 
conditions and relative momentum p = pi taken along the spin-quantization axis, i.e. the z-axis, and v Te \ = p//i is 
the magnitude of the relative velocity, p being the n-d reduced mass. The expression above is averaged over the spin 
projections m d = ±1, of the deuteron, however, the phase factor e mn = (— ) 1 / 2 ^ m « is ±i depending on whether the 
neutron has m n = ±1/2. 

The scattering states are calculated from realistic (strong interaction) Hamiltonians including two- and three- 
nucleon potentials by means of accurate hyperspherical- harmonics techniques @, 0, Hi , and are expanded in partial 
waves as 14] 

\pi;m n m d )^ = VI^J2 lL V2L+T(l/2m n ,lm d \ SJ Z ) (SJ z ,L0 \ JJ Z ) \ p; LS; J, J z )& , (2.2) 

LSJ 

where 

| p; LS; J, J z ) W = [1 T i J R(p)] t s ,L>s> I P;L'S';J,J Z ) . (2.3) 

L'S' 

In these equations, L denotes the relative orbital angular momentum between the two clusters, S their channel spin 
(either 1/2 or 3/2), and [ j R(p)]ls.l'S' is the (real) i?-matrix at center-of-mass energy p 2 /(2/i), from which phase- 
shifts and mixing angles in the channel specified by total angular momentum J are easily obtained. Lastly, the state 
| p; LS; J, J z ) is defined as in Eqs. (2.8) and (2.9) of Ref. [T3|, and in the asymptotic region (at large n-d separations) 
reduces to 



p;LS;J,J z ) -► -= ^2 w LS,L'S'&yi)[ Y L>(yi)®\iJk;S') 



ijk cyclic L' S 

w LS,L'S'{P^Vi) = $LL> 6 S s> jL>(pyi) + J RLs,L>s'(p)n L >(pyi) , (2.5) 

where | i, jk; S') denotes a state in which neutron i and a deuteron made up of nucleons jk are coupled to channel spin 
S' , and nL{x) are, respectively, the regular and irregular spherical Bessel functions, and = V{ — (r 3 + r^)/2 

is the n-d vector separation. 

We are interested in very low-energy neutrons, and thus it is sufficient to keep only S- and P-wave channels in the 
partial wave expansion of | pi; m n , mar • Equation (12. 1|) is then simply written as 



(2.4) 



-— = ^ V e m „ (1/2 m„, 1 m d \ JJ Z ) (1/2 m n , 1 m d \ SJ Z ) 

p dz ,/3v Iel ^ ^ 

(SJ Z ,10\ JJ z )1m[ ( -\p;lS;J, J z \ v pv \ p; OJ; J, J z ) (+) ] . (2.6) 
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The relation above follows by first noting that under time inversions, induced by the antiunitary operator T, the 
states in the convention adopted here transform as 

T \p;LS;J,J z )^ = (-) L+J -'> \p;LS;J,-J z )^ , (2.7) 

and then using the fact that the PV potential is i) a scalar under rotations (and hence its matrix elements are 
independent of J z ) and ii) invariant under time inversions, which lead to 

( -\p;0J;J,J z \v PV \ P ;1S;J,J Z ) ( ^ =- ( -\ P ;lS;J,J z \ v PV \ p; OJ; J, J z ) « . (2.8) 

In Eq. (|2.6[) the sums over S, J are restricted to S, J = 1/2 and 3/2, and thus the PV potential connects the doublet 
(quartet) S-wave state 2 Si/ 2 ( 4 S 3 / 2 ) to both the doublet and quartet P-wave states 2 Pi/2 ( 2 P 3 /2) an d 4 Pi/2 ( 4 P3/2)- 
Note that the additional contribution associated with the transition 4 S 3 / 2 — > 4 F 3 / 2 is neglected at the low energies 
of interest here. 



III. THE PARITY- VIOLATING POTENTIAL 



Two different models of the PV weak-interaction potentials are adopted in the calculations reported below. One 
is the model developed almost thirty years ago by Desplanques et al. 3] (and known as DDH): it is parameterized 
in terms of 7T-, p-, and w-meson exchanges, and involves in principle seven weak pion and vector-meson coupling 
constants to the nucleon. These were estimated within a quark model approach incorporating symmetry arguments 
and current algebra requirements [H, [Hj]. Due to the inherent limitations of such an analysis, however, the coupling 
constants determined in this way have rather wide ranges of allowed values. 

The other model for the PV potential considered in the present work is that formulated in 2005 by Zhu et al. Q 
within an effective-field-theory (EFT) approach in which only nucleon degrees of freedom are retained explicitly. At 
lowest order Q/A x , where Q is the small momentum entering the low-energy PV process and A x ~ 1 GeV is the scale 
of chiral symmetry breaking, it is parameterized by a set of twelve contact four-nucleon terms. In fact, it has been 
recently realized [161 ] that this set is redundant, and the number of independent contact terms can be reduced to five 
by Fierz-rearrangement. The derivation is summarized in Appendix |A"1 

The DDH and EFT PV two-nucleon potentials are conveniently written as 
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< °if > " = DDH or EFT 



(3.1) 
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where the parameters c" and operators O 



(n) 



r{n) 



are defined as 



1, . . . , 13, are listed in Table HI In this table the vector operators 



X 



(n) 



in) 



[Pij , fn{n 3 )] + , (3.2) 
X™. = i [ Pij , f n (rij)]_ , (3.3) 

where [...,.. .J^ denotes the commutator (— ) or anticommutator (+), and is the relative momentum operator, 
Pij = (pi — pj)/2. In the DDH model, the functions f x (r), x = tt, p and u>, are Yukawa functions, suitably modified 
by the inclusion of monopole form factors, 
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(3.4) 



In the EFT model, however, the short-distance behavior is described by a single function / M (r), which is itself taken 
as a Yukawa function with mass parameter p, 
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(3.5) 



with p ~ as appropriate in the present formulation, in which pion degrees of freedom are integrated out. 



In , the strong-interaction coupling constants of the 7T-, p-, and w-meson to the nucleon are denoted as 



while the weak-interaction ones as h\, h° hp, h'p 1 , hp, h^, where the superscripts 0,1, and 2 
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TABLE I: Components of the DDH and EFT models for the parity- violating potential. The vector operators X^L. and functions 



fx(r), x = n, p, ui, jj,, are defined in Eqs. (|3.2fl - (|3,3[l an< 4 Eqs. (|3.4[l ~ (|3.5fl . respectively. As outlined in the Appendix |A1 only 5 
operators and low-energy constants enter the pionless EFT interaction at the leading order, and in this paper they have been 
chosen to correspond to the rows 1, 4, 6, 8 and 9. . 



specify the isoscalar, isovector, and isotensor content of the corresponding interactions. In the EFT model, the five 
low-energy constants Ci,Cx, C2+C4, C5 and Cq completely characterize u|j- FT , to lowest order Q/A x (see Appendix |A"|) . 

In the analysis of the neutron spin-rotation observable to follow, we will report results for the coefficients I,^ DH and 
/ FFT in the expansion 



1# 

p dz 



13 

E 

71=1 



a ja 



(3.6) 



Thus we will not need to consider specific values, or rather range of values, for the strength parameters c". However, 
the 7" depend on the masses (and short-range cutoffs A x for the DDH model) occurring in the Yukawa functions. 
Two different sets of values are used for A^ and fi in the present study, as listed in Table HT1 Both were used in the 
DDH-based calculations of PV two-nucleon observables in Refs. P, 0). 



IV. CALCULATION 



The relevant matrix elements which need to be calculated in Eq. 



are given by 



p; IS; J, J z I v PV I p; 0J; J, J,)W = 3 / dxdy ^ 8 > JJt (x, y) 



(+) _ 



.(-) 



(+) 



23 Y0J;JJ : 



( x .y) , 



(4.1) 



where x = r2 — r,3 and y = ri — (r2 + r^)/2 is the n-d vector separation defined earlier. The factor of 3 on the 
right-hand-side simply accounts for the three identical contributions associated with the sum over pairs included in 





A,r 


Ap 


A^ 




A 4 


DDH-I 


1.72 


1.31 


1.50 


EFT-I 


0.138 


DDH-II 


00 


00 


00 


EFT-II 


1.0 



TABLE II: The two different sets of values, in GeV units, for the short-range cutoffs A x (mass fi) used in the DDH (EFT) 
model. Note that the masses m n , m p , and ra w are taken respectively as 0.138, 0.771, and 0.783 in units of GeV. The row with 
A x — 00 corresponds to point-like couplings. 
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PV 

v , 



PV 

V 



where a=DDH or EFT. 

The (antisymmetric) wave functions i/»(x, y), dropping subscripts and superscripts for brevity, are written as 

^(x,y)= 22 i>i;jk(*i,yi) , (4.3) 

ijk cyclic 

where 

Xi = Yj -r k , yi = ri - (rj + r k )/2 , (4.4) 

(ijk)=(123), (231), and (312) are the three cyclic permutations of the particles, and Xj = x and yj = y. The 
decomposition given in Eq. (|4. 3|) corresponds to the three possible partitions of three nucleons in 1+2 clusters — 
indeed, the asymptotic wave functions in Eq. (|2.4|) have precisely this form. Note that in the integral of Eq. (|4.1|) 
there are terms like (again in a schematic notation) 

J dxdy^t. 23 (xi,yi)«23^i;23(xi,yi) , (4.5) 

corresponding to the partition nucleon 1 and nucleons 23 in the initial and final n-d states, and with the PV interaction 
acting on nucleons 23. These n-d continuum wave functions have "core" parts which vanish at oo, and asymptotic 
parts oscillating in the y\ variable. The integrals involving core parts of either the left or right (or both) wave functions 
are converging, while the integral involving the left and right asymptotic wave functions vanishes, since the integrand 
is odd under inversion yi — > — yi. The remaining partitions give converging contributions to the integral, since the 
deuteron cluster in the left and right wave functions is made up of different nucleons and V'2;3i( x 2) y2)j V'3;12( x 3j Yz) ~ > 
as yi — > oo. 

The wave functions for an assigned spatial configuration (x, y) are expanded on a basis of 8 x 3 spin-isospin states 
for the three nucleons as 

24 

^(x,y)=5^Va(x,y)|o), (4.6) 

a=l 

where the components ip a (x,y) are generally complex functions, and the basis states | a)= | (n i)i{p i)2{n !) 3 ), 
I (ni)i(n l)2(j>i)3), and so on. Matrix elements of the PV potential components are written schematically as 

24 „ 

(f\0\i)= / dxdy V$,„(x,y) [0(x)] a6 ^, b (x, y) , (4.7) 

a,b=l J 

where [0(x)] afc denotes the matrix representing in configuration space any of the components in Table fl] Note that 
the operators X23 \r occurring in are conveniently expressed as 

X23?- = , Xg + = -i [2 f n {x) V x +x/;(x)] , (4.8) 

where the gradient operator V K = (V2 — Va)/2 acts on the right (initial) wave function, and f'(x) = df(x)/dx. 
Gradients are discretized as 

ip(...r i +de a ...)-i/)(...ri-6e a ...) 
Vi, a V(x,y) ^ y§ ' ( 4 - 9 ) 

where S is a small increment and e Q is a unit vector in the a-direction. Matrix multiplications in the spin-isospin space 
are performed exactly with the techniques developed in Ref. [TtJ • The problem is then reduced to the evaluation of 
the spatial integrals, which is efficiently carried out by a combination of Monte Carlo (MC) and standard quadratures 
techniques. We write 

r 1 N < 

{f\0\i)= dxdy^(x,y)^— ]TF(c) , (4.10) 

J c c=l 
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where the c's denote uniformly sampled directions (x, y) (total number jV c ). For each such configuration c, the 
function F is obtained by Gaussian integrations over the x and y variables, i.e. 

F(c) = (4^) 2 ]T / dxx 2 d 2/2/ 2 ^. Q (x, y )[0(x)] Qb ^ b (x,y) . (4.11) 

Convergence in the x and y integrations requires of the order of 40 Gaussian points, distributed over a non-uniform 
grid extending beyond 40 fm, while N c of the order of a few thousands is sufficient to reduce the statistical error in 
the MC integration at the percent level. Thus, the present method turns out to be computationally intensive. In 
this respect, we note that a direct MC evaluation of the six-dimensional integral in Eq. (|4.7|) by Metropolis sampling 
from a probability density function W(x, y) cx l^o^y) | 2 > where ipo is the triton bound-state wave function, turned 
out to be impractical. This was also the case for other choices of W, such as W(x,y) oc x x y" | 0(x) | 2 | 4>(y) | 2 , where 
4> is the deuteron wave function. The reason for this difficulty can be appreciated by examining Fig. [TJ in which the 
integrand associated with the pion-range component of the DDH potential, specifically the matrix ele ment of 
between the 2 Si/2 and 2 Pi/2 n-d states at zero energy, is plotted as function of the hyper-radius p = ^/ x 2 + Ay 2 /3. 
Note the change of sign in this function, due to the node in the continuum 2 Si/2 state (thus ensuring its orthogonality 
to the 3 H bound state), and its long-range character. This last feature, in particular, makes it difficult for an efficient 
sampling of the large p region with the importance functions W mentioned above. 

The computer codes were successfully tested by carrying out a calculation based on Gaussian wave functions for 
the initial and final states, as described in the following subsection. 



A. Test calculation 



We have performed calculations with Gaussian-type wave functions, for which the matrix elements can be computed 
analytically. In Eq. (|4.ip . we have replaced the realistic n-d wave functions with the following ones 

^LS-JJ. (X, y) = ^-JL^ n LS ;JJ z , (4.12) 



/4tt 

with 



n LS ;j, h = 22 Vi F LsAVi) [Y L (yi) ® [siSjkjsljj [tiT jk ] TTz , (4.13) 

ijk cyclic 



where p is the hyper-radius (independent of the particle permutation considered), Si (U) denotes the spin (isospin) 
state of particle i, and Sjk (Tjk) the spin (isospin) state of particles jk. We have considered only states of total 
isospin T, T 2 =l/2,-l/2, appropriate to describe the n-d channel — possible admixtures of T=3/2 components induced 
by isospin symmetry breaking terms in the strong and electromagnetic interactions are expected to give negligible 
contributions to PV observables. 

The quantum numbers Sjk and Tjk are respectively taken to be 1 and 0, while the factors F^sjiyi) selected for the 
various channels are reported in Table lilTl Note that FLs,j{yi)—^ except for L—0 and S=J—3/2, since in this case 
the spin part [siSjfc] 3 y 2 is totally symmetric, and it is impossible to construct a fully antisymmetric wave function 

with only the isospin part. Hence, the choice F(yi)~y 2 . In all other channels, non-vanishing functions flLS,JJ z & r e 
obtained by simply taking i 7 LSj(?/i) = l- 

The analytical calculation is simplified considerably by expressing the function Qls-,jj z i which involves the three 
set of coordinates (xi,yi), (x2,y2), and (x3,ya) corresponding to the three partitions (1,23), (2,31), and (3,12), in 



S 



j F LSJ { yi ) 
~t 1 



TABLE III: Quantum numbers and factor Flsj(ih) entering the Gaussian wave functions given in Eq. (14.121) . 



10 20 30 40 50 

p[fm] 



FIG. 1: The function G(p), defined so that dpG(p) is the matrix element of the operator O' 1 ' between the 2 Si/2 and 2 Pi/2 
n-d states at zero energy, plotted as function of the hyper-radius p = x 2 + 4 y 2 /3. 

terms of linear combinations of functions of, say, the first set of spatial variables, i.e. x = xi and y = yi as in 
Eq. (14. ip , and spin-isospin states of pair 23 only. We find: 

^LS;JJ Z = c nf,n„,^,S x ,j x ,i y ,j v ,T x ^n x ,k y / x ,S x ,3 x ,iy,3v,T x ' ( 4 ' 14 ) 

Tlx -i^ly , S x ijx i^-y i jy j Tx 



where 



r>JJ z ;TT z _ l x +n x „.i y +n v 



[Y^)S x ] jm ® [^(y) Sl ]J ^ [T x h] TTz , (4.15) 



and n x are n y are non-negative integers, S x and T x are the spin and isospin states of pair 23, and c^ x J ny £ x Sx Jx ly j Tx 
are numerical factors, which can easily be expressed in terms of Wigner coefficients. 

It is now relatively simple to compute matrix elements of the operators in Table [J between states 
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fl n ^ i s j i j t ^ making use of the following relations 



J die [Y ex (x)S x ]l mx (a 2 + a 3 ) ■ x [Y,, (^J.^ = -2 J 3x+ ^^' 2 6t.,i. , (4.16) 

and similar ones for the spin-space and isospin operators occurring in the potential components 2 3^ • ^ n particular, 
those for (<r 2 — <x 3 ) • and (er 2 x cr 3 ) • x are listed in Sec. III.E of Ref. The remaining task of integrating over 
the radial variables x and j/ is carried out accurately by means of Gauss quadrature formulae. 

The results of the "analytical" calculation just described are found to be in agreement, for each of the 13 components 
of the PV potential, with those produced by a MC calculation using the same input wave functions as in Eq. (|4.12|) . The 
spin-isospin algebra is performed as outlined in the previous section; however, because of the short-range character of 
the wave functions in this case, spatial configurations turn out to be efficiently sampled with the Metropolis algorithm 
by a probability density proportional to exp(— 2 (3p 2 ). In a typical calculation, the values (3 = 0.10-0.25 fm -2 have 
been selected. 



V. RESULTS AND CONCLUSIONS 



The results for the coefficients I" in Eq. (|3.6|) . obtained with n-d continuum wave functions corresponding to the 
AV18 and AV18/UIX strong-interaction Hamiltonians, are reported for the DDH and pionless EFT PV potentials in 
Tables IIVI and [Vj respectively. The labels I and II for the columns in each of these tables refer to the two sets of 
cutoff parameters, as given in Table HT1 The subscript n in 1° specifies the operators as listed in Table |U Note that 
results for n—li are not re por ted, since the coupling constant hp' multiplying these contributions has been estimated 
to be considerably smaller [15| than the "best" values estimated for the other coupling constants in the original DDH 
reference Q. 





DDH-I 


DDH-II 


n 


AV18 


AV18/UIX 


AV18 


AV18/UIX 


1 


0.256E+03 


0.270E+03 


0.257E+03 


0.274E+03 


2 


-0.444E+01 


-0.691E+01 


-0.719E+01 


-0.118E+02 


3 


0.444E+01 


0.401E+01 


0.732E+01 


0.761E+01 


4 


-0.231E+01 


-0.881E+00 


-0.332E+01 


-0.148E+01 


5 


-0.247E+01 


-0.122E+01 


-0.387E+01 


-0.234E+01 


8 


0.420E+01 


0.362E+01 


0.543E+01 


0.516E+01 


9 


0.111E+01 


-0.117E+00 


0.136E+01 


-0.189E+00 


10 


-0.253E+01 


-0.991E+00 


-0.314E+01 


-0.141E+01 


11 


-0.280E+01 


-0.144E+01 


-0.369E+01 


-0.225E+01 


12 


0.316E+01 


0.339E+01 


0.455E+01 


0.546E+01 



TABLE IV: The coefficients 7° DH (in fm) in Eq. (|3~6)l corresponding to the DDH PV potential in combination with the AV18 
and AV18/UIX strong-interaction potentials for the two different sets of cutoff parameters given in Table [II] The statistical 
errors associated with the Monte Carlo integrations are not shown, but are typically at the 1-2% level. Note that the coefficients 
Z^?6 7 vanish because of isospin selection rules. 

A quick glance at Table HVl makes it clear that the matrix element of the long-range component of the DDH potential 
due to pion exchange is almost two orders of magnitude larger than that of its short-range components induced by 
vector-meson exchanges. It is also clear that this contribution is fairly insensitive to the choice of cutoffs as well as 
of input strong-interaction Hamiltonians used to generate the n-d S- and P-wave channel states. In particular, one 
should note that the AV18 and AV18/UIX models predict rather different values for the scattering length in the 2 Si/ 2 
channel [8]. On the other hand, predictions for the 4 Si/2 scattering length are very close to each other. This is reflected 
in the contributions to /p DH due to the doublet and quartet S-wavc channels. Wc find that for the case of DDH-I, for 
example, these contributions for the AV18 and AV18/UIX are, respectively, 45.0 fm and 74.2 fm in the 2 S!/ 2 channel, 
and 210.6 fm and 196.6 fm in the S j/g channel, which add up (up to truncation errors) to the values 256 fm and 270 
fm, reported in the first row of Table ll VI We observe that the doublet scattering length calculated with the AV18/UIX 
is very close to the experimental determination of this quantity, while the quartet scattering length, which is much 
less sensitive to the presence of three-nucleon interactions, is found to be in agreement with the experimental value 
with both the AV18 and AV18/UIX models. Hence, it is reasonable to expect that any strong-interaction Hamiltonian 
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that reproduces the experimental scattering lengths should lead to predictions for the long-range component of the 
DDH potential, which are very similar to those reported here for the AV18/UIX model. 

There is considerable model dependence in the results obtained for the individual contributions due to vector-meson 
exchanges. However, this model dependence has little impact on the neutron spin rotation, as it can be surmised 
from Table IVI1 where we list predictions obtained for this observable. We adopt two different sets of strong- and 
weak-interaction coupling constants, namely those reported in Tables I and II of Ref. Q . The set of coupling constants 
denoted as DDH-best (from Table II of Ref. @) is the "best value" set as given in the original work . Recently, the p- 
and cj-meson weak-coupling constants were adjusted in Ref. [l[ to reproduce precise measurements of the longitudinal 
asymmetry in p—p elastic scattering (see Table I of Ref. Q). The resulting set of coupling constants is denoted as DDH- 
adj. Note that both sets use /i* =4.56 x 1CP 7 . Assuming a liquid deuterium density of p = 0.4 x 10 23 atoms cm -3 , we 
obtain d0/dz » 0.53 x 10" 7 rad-cm" 1 (0.56 x 10~ 7 rad-cirr 1 ) with the AV18/UIX+DDH-best (AV18/UIX+DDH-adj) 
model. About 90% of this prediction comes from the long-range component of DDH, while its short-range components 
contribute the remaining 10%. In particular, the reason for the larger short-range contribution obtained with DDH- 
adj relative to DDH-best lies in the fact that this contribution is dominated by the n — 3 term in Eq. (|3.6j) , for which 
the combination of coupling constants is C3 = —g p h° p (1 + K p )/m (see Table Jj). In DDH-adj C3 is roughly a factor of 
two larger than in DDH-best because of different values for h° p and the tensor coupling constant n p . 

It is interesting to note that the neutron spin rotation in ri-p scattering is predicted to be about one order of 
magnitude smaller than in the present case |2j, which should make its measurement in n-d scattering considerably 
easier, at least in principle. 



EFT-I EFT-II 



n 


AV18 


AV18/UIX 


AV18 


AV18/UIX 


1 


0.257E+03 


0.274E+03 


0.392E+01 


0.706E+01 


4 


-0.154E+03 


-0.508E+02 


-0.127E+01 


-0.682E+00 


8 


0.260E+03 


0.189E+03 


0.233E+01 


0.251E+01 


9 


0.244E+00 


-0.359E+02 


0.557E+00 


-0.955E-01 



TABLE V: Same as in Table IIVI but for the pionless EFT PV potential. Note that there are no potential components with 
Oy with n—2, 3, 5, 7, 10 and 11, and that the coefficient In=e vanishes because of isospin selection rules. 

The coefficients J„ for the relevant operators entering the pionless EFT PV potential, that is n=l, 4, 6, 8, and 
9, are reported in Table fVl Note that the n=6 operator is isotensor, and therefore its matrix element vanishes here. 
The /^ft coefficients are all of the same order of magnitude, since the radial functions are taken to be the same for 
all n, fn FT (r) = /^(r). Of course, they depend significantly on the value of the mass p — either p = m, (EFT-I), as 
appropriate in the present pionless EFT formulation, or p — 1 GeV (EFT-II), the scale of chiral symmetry breaking, 
as appropriate in the formulation in which pion degrees of freedom are explicitly retained. Indeed, in this latter 
formulation the leading order component of u PV has the same form as the pion-exchange term in DDH. 

We note that the coefficient I^ZJ 1S weakly dependent on the choice of input strong-interaction Hamiltonian — AV18 
or AV18/UIX — in the case of the EFT-I model. On the other hand, this is not the case the for I^EIs 9 coefficients. 

Finally, rough estimates have been made for the range of values allowed for the low-energy constants C\, C2 + C4, 
C5, Ci, and Cq in Ref. [5j. However, at the present time a systematic program for their determination is yet to be 
carried out. In view of this, we refrain here from making EFT-based predictions for the spin rotation observable. 



AV18 AV18/UIX 





DDH-best 


DDH-adj 


DDH-best 


DDH-adj 


7T 


0.457 


0.457 


0.485 


0.485 


p-UJ 


0.063 


0.092 


0.046 


0.075 


TOT 


0.520 


0.550 


0.531 


0.560 



TABLE VI: Spin rotation for ri — d scattering in units of 10 -7 rad-cm _1 at zero energy, obtained for the DDH-I model assuming 
a liquid deuterium density of p = 0.4 x 10 23 atoms cm -3 . The columns labeled DDH-best (DDH-adj) list the spin rotation 
predicted using the "best" ("adjusted") values for the weak-interaction coupling constants, as reported in Table II (I) of Ref. 0. 
The row labeled 'V ("p-w") lists the results obtained by including only the pion (p and uj mesons), while that labeled "TOT" 
gives the total contributions. 

In conclusion, we have calculated the neutron spin rotation, induced by PV components in the nucleon-nucleon 
potential, in n-d scattering at zero energy. We find that this observable is dominated by the contribution of the long- 
range part of the PV potential due to pion exchange, and that the model-dependence of this contribution is weak. The 
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predicted d(j)/dz is ~ 0.5 x 10~ 7 rad-cm -1 — an order of magnitude larger than expected in n-p scattering 0. Thus a 
measurement of the spin rotation in n-d scattering could provide a further constraint, complementary to that coming 
from measurements of the photon asymmetry in n-p radiative capture, on the strength of the long-range component 
of the hadronic weak interaction. 
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APPENDIX A 

The leading order PV contact Lagrangian has been written in Ref. [a] as consisting of twelve operators, 

6 

£ = ^(0*0* + 0A) , (Al) 
i=i 

where 

01 = '(/^'V^/asV' , 0} = ^17^75^ (rpla^ip) , 

2 = -017 / VV'T37m75V' , 2 = i>\^^i)d v (ipT 3 a^'iJj) , 

3 = f T a r^T a ^ l5 ^ , o 3 = ^^v^tv^) , (A2) 

Oi = ^37^-017^75^ , 4 = ^T 3 j^-f 5 ipd u (ipla^ip) , 

6 = ie ab3 i{;T a ^T b j 5 ^ , O e = ie afc3 Vn- a 7^^(^T% i ,7 5 V') , 

where the nucleon field ip is a doublet in isospin space, and the operators 1 and t° above act on this space. In addition, 
Tab is a 3 x 3 diagonal matrix with elements (1^1,-2). Using Fierz transformations and fields' equations of motion, 
the following relations may be shown to hold [lq . 

O3 = Ox , 

2 - 4 = 20 6 , 

6 3 = 2m(0 1 + 3 )-30 1 , 

6 2 + Oi = m(0 2 + Oi) , [A6) 
2 -Oi = -2m0 6 - 6 , 
65 = 5 , 

so that the number of independent operators can be reduced to six. Moreover, since the relativistic operators Oq 
and 6 give rise to the same structure in the leading order of the non-relativistic expansion, as already observed in 
Ref. [5j, the PV two- nucleon non-relativistic contact Lagrangian consists of only five operators, 

C PV:NN = —{c l (N^aN-NhvN-N^NNhv-vN)~C l e l:jk N 1 i<T l NV : >(N 1 i<T k N) 

x 

-(0 2 + Ci)t ljk [N* t 3 a l NV 3 (N^ a k N) + A* a 1 AV^A* t 3 a k A)] 

~C 5 l ab e ljk N^r a a i NV j (N^ r b a k N) +0 6 e ab3 V(AV A) ■ AVoW} , (A4) 

where the notations of Ref. [f| have been chosen for the coupling constants. Accordingly, the pionless EFT potential 
in coordinate space reads as 

eft _ ^[i 2 

ij A3 



01 + (0 2 + 4 ) ( -^"^ ) + 05 l ab T? Tj 



+1 0! fa x CTj) ■ [ P y , f,(r)}_ + i C 6 e ab3 T«T b (cr z + a,) ■ [ Pij , / M (r)]_} , 
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where Py is the relative momentum operator and /^(r), defined in Eq. (|3.5p . is introduced to regularize the potential 
at short-range. 
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